Phase-dependent quasiparticle tunneling in Josephson junctions: Measuring the 
cos(/?-term with a superconducting charge qubit 
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We investigate quasiparticle tunneling in a Cooper-pair box which is embedded in a supercon- 
ducting ring to allow control of the total phase difference across the island. The phase affects 
the transition rate between different electron number parity states of the island, which can be ob- 
served in experiment by established means. The phase dependence also leads to what is known 
as the cos y>-term in the tunneling characteristics of classical Josephson junctions. This effect has 
remained controversial for decades; the proposed scheme opens an independent way to probe it. 
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Introduction. Since the discovery of the Joseph- 
son effect 1 ' 2 it has been known that, in addition 
to the phase-dependent Cooper-pair current and the 
well-known quasiparticle current, there exists a phase- 
dependent quasiparticle- pair interference term 1 ''. It is 
often called the cos <^-term since in the simplest approxi- 
mation it adds a dissipative term proportional to cos tp to 
the RCSJ model of a current-biased Josephson junction, 

1 2e 2e 

Ctp + — (1 + ecostp)tp + — Ij sirup = — I. (1) 
R n h 

Here C is the junction capacitance, R the subgap re- 
sistance, Ij the critical current, and / the bias current. 
The theory, based on the BCS and tunneling Hamilto- 
nian predicts for slow variations of tp the value e ~ 1 , 
whereas various experiments suggested e ss — 1 1-8 . How- 
ever, these experiments were difficult to interpret, and it 
is probably fair to say that the discrepancy and its ori- 
gin was never fully resolved. Recently, after a period of 
reduced interest, the role and properties of quasiparti- 
cle tunneling in Josephson junctions has again attracted 
attention, since it provides a possibly important contri- 
bution to the relaxation in superconducting qubits . 

In this article we show that the quasiparticle-pair inter- 
ference term also affects the single-electron tunneling rate 
in low-capacitance superconducting tunnel junctions. We 
propose a scheme to measure it via transitions between 
different charge eigenstates of a Cooper-pair box (CPB) 
embedded in a superconducting ring as shown in Fig. la. 
The system is particularly sensitive to quasiparticle tun- 
neling, as these processes change the particle-number 
parity of the island. The signature of the quasiparticle- 
pair interference is a transition rate with additional terms 
proportional to e cos tpi , where the phases tpi can be con- 
trolled by the applied magnetic flux, allowing tuning the 
effect 11 . States with even or odd parity can be dis- 
tinguished by measuring the effective quantum capaci- 
tance of the CPB, as demonstrated in the experiments of 
Ref. 13. This technique has already been used to probe 
quasiparticle tunneling in the presence of their typical 
non-equilibrium density, but no attention has been paid 
to the phase dependence and the interference term 14 . 




FIG. 1: (a) A Cooper-pair box between two Josephson junc- 
tions in a superonducting loop threaded by a flux <& x . The 
small island is capacitively coupled to a voltage gate U . (b) 
The quasiparticle current 7 qp and the quasiparticle-pair inter- 
ference term 7j2 for a voltage-biased Josephson junction with 
tunneling resistance _Rt at T = 0.9T C . 



One of the difficulties of the experiments aimed at re- 
solving the cos <^>-term lies in the fact that with current- 
biased junctions one can control and probe the phase- 
dependent quality factor of small oscillations only for 
a restricted range of phase values. SQUID configura- 
tions ' ' where the flux can be biased and the conductance 
can be measured over the whole phase periodicity offer, 
in principle, a solution to this problem. However, such se- 
tups suffer from the problem that variations in tp change, 
simultanously with the dissipation due to the cos c/3-term, 
also the oscillation frequency. This in turn can alter the 
effect of the electromagnetic environment, which in many 
cases is the main source of dissipation 1 . Our proposed 
measurement scheme overcomes these problems as for the 
CPB the dominant energy scale, the charging energy, is 
independent of the flux, and for the change of particle- 
number parity of the island no significant process com- 
peting with quasiparticle tunneling exists. 

A possible reason for the sign discrepancy of e lies 
in the dependence of the quasiparticle-pair interference 
term on the voltage V cx tp. For constant voltage it can 
be easily evaluated, the result 7j2 is plotted in Fig. lb. It 
changes sign at a voltage matching the superconducting 
gap, eV = 2A. Accordingly, e is expected to be positive 
below the gap and and negative above. A broadening of 
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FIG. 2: (a) A quasiparticle tunneling across the Josephson 
junction from left to right gains the even-odd energy differ- 
ence SE. It can involve either electron tunneling from left to 
right (upper process), or electron tunneling from right to left 
(lower). In the latter process a Cooper pair breaks on the 
right and the tunneling electron recombines with an unpaired 
electron on the left, (b) Breaking of a Cooper pair during 
quasiparticle tunneling. This process sets in when SE > 2A 
and has also two possible electron tunneling directions. 



the quasiparticle states can lead to a negative sign also 
at subgap voltages 10 '' 1 ' . However, for typical materials 
used today this effect should be negligible. In the setup 
of Fig. la, without further modification, only thermally 
excited or non-equilibrium quasiparticles can tunnel (see 
Fig. 2a). This is a subgap process corresponding to pos- 
itive sign for e. In order to probe the sign change of e 
we propose in the second part of this article a modified 
setup with a lower superconducting energy gap in part 
of the ring. Parameters are chosen such that enough en- 
ergy is available to split a Cooper-pair (see Fig. 2b) in 
a quasiparticle tunneling process across one of the JJs. 
This corresponds to a process above the gap and, hence, 
to a negative sign for e. 

System. The system considered is shown in Fig. la. It 
consists of a small island (Cooper-pair box) capacitively 
coupled to a gate and embedded via two Josephson junc- 
tions in a superconducting ring that is threaded by an 
external magnetic flux § x = (H/2e)<fi x . If the inductance 
of the ring L is small, we can assume that the phase dif- 
ference across the split CPB is equal to the phase bias 
4> x . The coherent time evolution of the system is then 
described by the Hamiltonian 

ffcoh = Ec(N-N G ) 2 -E n cos(<p+0 x )-Ej 2 cos(<p), (2) 

where Eji are the Josephson couplings of the two junc- 
tions, Ec = e 2 /2Cs denotes the charging energy of the 
island with Cs = C\ + C 2 + Cq being the sum of the ca- 
pacitances surrounding the island, and Nq = —CoJJje. 
The phase ip and the number of excess electron charges 
on the island N satisfy the periodic commutation rela- 
tion [N/2,e ±lip } = ±e ±l<p . Cooper pairs tunnel coher- 
ently across the two junctions with a relative phase shift 
4> x , the phase difference across the split CPB. It can be 
treated classical variable. 

In the following we assume that Ec > Eji,Ej 2 , and 
the gate is biased at Nq = 1. Therefore we restrict our 
analysis to the subspace spanned by the island charge 
states |0), |1), and |2), with N = 0, 1, or 2 extra electron 
charges. Choosing the energy of the state |1) as reference 



the coherent part of the Hamiltonian takes the form 
#coh = £c(|0)<0| + |2)<2|) 

- (e^ rff |2)(0| + e-^ ff |0)(2|) . (3) 
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Ej2 , with 



Here w e introduced Eje llPaii 
Ej = yjE 2 n + £j 2 + 2E n E J2 cos <f> x > and tan^ cff = 
sin X / (cos 0a +Ej 2 /E n ). If E, n > E J2 the effective 
flux ifeff spans the whole range from — n to n. The Hamil- 
tonian (3) does not connect states with even and odd 
numbers of electrons on the island. It has the eigenstates 
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I) = ^=(|0>+e^|2», 



(4) 



in addition to the "odd" state |o) = |1), with eigenen- 
ergies E t = E c + E 3 /2, E i = E c - E 3 /2 and E Q = 0. 
The relevant energy difference in the following is the one 
between the two lowest eigenstates SE = E^. 

The quasiparticles are described by the BCS Hamil- 
tonian. Their tunneling across the junctions leads to 
transitions between states having even or odd number 
of electrons on the island. For junction i the tunneling is 
described by the Hamiltonian 

H Tl = U ^[(w k 7kt + v k7k4.)(wi7i t + v\Tn)Ti 



kl 



(5) 



where 7^ is a quasiparticle annihilation (creation) oper- 
ator of the state 1 and spin a in the island. The states in 
the loop side of the junctions are labelled as k. The 
quasiparticle operators are related to the correspond- 
ing electron operators via the Bogoliubov transformation 

4f = Wi7it + and -U = ~ v nit + M i7iV Hcrc thc 
BCS coherence factors u and v are real numbers as the 
phase dependence is included with the charge-transfer 
operators f 2 = \2){l\ + |1)(0| and T\ = e l ^/ 2 f 2 , which 
take care of the corresponding island-charge changes in 
the tunneling events. They have similar (but halved) 
external flux dependence as the Cooper-pair tunneling, 
described by Eqs. (2-3). The average tunneling matrix 
elements U are related to the normal state tunnel resis- 
tances via i?Ti = h/t 2 D 2 ire 2 , where D is the density of 
states at the Fermi surface, including spin. 

Tunneling processes. The tunneling Hamiltonian 
(5), when treated perturbatively in leading order, in- 
cludes three distinct types of quasiparticle tunneling pro- 
cesses. In the first a quasiparticle, either thermally ex- 
cited or a non-equilibrium one, tunnels across the junc- 
tion. The process can still involve two electron-tunneling 
directions (terms like 7 k 7iT or 7 k 7iTt), as visualized in 
Fig. 2a. The second process needs at least an energy 2A 
to break a Cooper pair, which then provides an electron 
that tunnels across the barrier (terms like 7^7/ T and 
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7 k 7i"Tt illustrated in Fig. 2b). For the moment we can 
ignore it, but it plays a role in the modified design dis- 
cussed later. The third tunneling type corresponds to the 
case when an electron tunnels across a junction, with a 
simultaneous recombination of two excitations originat- 
ing from different sides of the junction (terms like 7k 71 T 
and 7k7il^)- The process does not contribute under the 
operating conditions considered here. 

The even eigenstates are superpositions of the two 
charge states N = and iV = 2, as given by Eq. (4). 
Therefore, in the transition from, e.g., | \.) to |o) both 
electron tunneling directions contribute. They are coher- 
ent, and hence the two amplitudes have to be added be- 
fore squaring in the Golden rule calculation. This intro- 
duces an extra contribution to the transition rate, which 
is proportional to the cosine of the relative phase. It cor- 
responds to the cos yj-term in single Joscphson junctions 
mentioned in the introduction. This contribution would 
not be contained in the standard semiconductor model 
treatment 15 . Moreover, the interference term depends 
on which tunneling process is dominant. It differs for the 
tunneling of an existing quasiparticle from the process 
which involves breaking of a Cooper-pair, the difference 
being in the sign of e of the interference term. 

Transition rates. In the following we assume Ej ^> 
fceT, which implies that of the two even parity states 
only the lower eigenstate | \.) needs to be considered. The 
quantity to be calculated is the transition rate Tt from 
this state to the odd ground state |o). The opposite tran- 
sition from odd to even occurs either through tunneling of 
thermally excited quasiparticles or of a non-thermalized 
quasiparticle, as discussed in Ref. 9. These two odd-to- 
even processes also show interference, but the rate de- 
pends strongly on the energy barrier 8E, which can have 
a further dependence on the magnetic flux through Ej. 
In order to avoid this complication and to observe only 
the interference effect we concentrate on the even-to-odd 
transition. It can be resolved in the experiment where the 
transitions are seen in real time as jumps in the phase of a 
reflected weak microwave signal applied to the gate, due 
to simultaneous changes in the quantum capacitance of 
the system 14 . One of the first quantum capacitance ob- 
servations has been achieved using the exact same setup 
we consider in this paper 21 . 

The total even-to-odd transition rate can be decom- 
posed into four contributions, 
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where r^ 2 ^ is the forward/backward rate for quasiparti- 
cle tunneling through junction 1 and 2, respectively. Ex- 
panding the time evolution of the density matrix up to 
second order in the tunneling Hamiltonian (5) we obtain 
the forward tunneling rate for junction 1 
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is the quasiparticle excitation 

energy and /l(i) [£&(;)] the excitation occupation probabil- 
ity in the loop (island) side of junction 1. Similar expres- 
sions apply for the backward tunneling rate and junction 
2. The information of the interference is contained in the 
function 
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where ipi = (f> x — y ff and Lp 2 = <p e s. The first term 
corresponds to the quasiparticle current, while the coscp- 
term is the quasiparticle-pair interference : . The last term 
gives no contribution if the distribution is symmetric with 
respect to electron (k > kp) and hole (k < kp) branches, 
which is the case for thermal distribution. 

To proceed we use a quasiparticle distribution 
fsn(E) = f e q(E — Sfi) where 5 fx is a shift in the chemi- 
cal potential and f cq (E) the equilibrium Fermi function. 
This description in terms of a shifted chemical potential 
implies that all quasiparticles arc in thermal equilibrium, 
but their density is out of equilibirum (5[i > 0). This 
type of non-equilibrium distribution may emerge, e.g., in 
the presence of an electromagnetic radiation 1 *, and it has 
been realized in recent experiments 14 ' 19 . It has also been 
used to model quasiparticles in the odd parity state 20 . 
We further assume that the quasiparticles are distributed 
uniformly in the system and, hence, 

temperatures the decay rate (6) then reduces to 
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where r = 5E/A, e = 1/(1 + r), and n qp = 
2D dEfs^E 1 'v 'E 2 — A 2 is the quasiparticle density. 

The transition rate Tt as function of the applied 
phase <j) x for typical experimental conditions is plotted 
in Fig. 3a. We have chosen an asymmetric situation with 
mm{E } {$ x )}/k B = (E. n - E J2 )/k B » 0.2 K, and the 
condition Ej 3> k-gT is satisfied for the chosen tempera- 
ture T = 50 mK. We observe an almost sinusoidal varia- 
tion of the total rate, which arises from the interference 
effect. From the individual contributions of the junc- 
tions we see that this pattern originates from tunneling 
across the weaker junction 2. This is the case because for 
Ej\ 3> Ej2 one has ip e s ~ 4> x , and for quasiparticle tun- 
neling across the junction 1 the first two terms in Eq. (8) 
always show destructive interference. For junction 2 the 
interference is proportional to ~ 1 — e cos <p x , resulting in 
the sinusoidal phase dependence. We emphasize that the 
phase dependence of the total rate originates predomi- 
nantly from the cos </?-terms (quasiparticle-pair interfer- 
ence). A competing effect could be the weak flux depen- 
dence of the energy level 5E = Eq — Ej/2 through Ej. 
However, this contribution remains always small because 
of the high charging energy of the CPB. The absence of 
the interference would correspond to e = 0. 

The temperature dependence of the rates is governed 
by the quasiparticle density n qp . At high temperatures 
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FIG. 3: (a) The transition rate Tt for the even-to-odd par- 
ticle number switching as a function of the phase bias <j> x . 
Also plotted are the individual contributions for junction 1 
(ri = + rj") and junction 2 (T2). At the considered low 
temperature, T — 50 mK, the transition is dominated by non- 
equilibrium quasiparticles, for which we assume the density 
?iqp = 10/^rn 3 , similar to recent experiments 14 . We use typi- 
cal aluminum density of states DA = 2.8 x 10 6 /^im 3 , and as- 
sume the Ambegaokar-Baratoff relation between Eji and R^i . 
Other parameters are En = 4E J2 = A/10 = E c /5 = 20 /xeV. 
(b) The rates for the modified system with a reduced energy 
gap in parts of the loop. The tunneling rate in junction 1 does 
not depend on the applied flux and the rate of junction 2 has a 
sinusoidal variations, but with opposite oscillation phase and 
hence sign of e as compared to (a). The parameters are the 
same as in (a) except for Ej2 = 5 * 10~ J /^eV, Ec = 150 fJ.eV, 
and Ai = 4A 2 = A. 

it is proportional to exp(— A/fcgT), at low temperatures 
it saturates due the assumed shifted chemical potential 
<5/i > 0. Well below this cross-over the individual junc- 
tion rates grow (weakly) with decreasing temperature for 
a constructive interference and decrease for a destructive 
one. However, for a system with clearly asymmetric junc- 
tions the total rate always increases with temperature, 
since for the larger of the junctions the interference stays 
practically destructive for all <p x . 

Sign change of the interference term. We now 
extend our discussion to situations where the change 
of sign in e can be probed. For this purpose we as- 
sume that the loop consists of two different supercon- 
ducting materials, one having a much lower energy gap 
than the other. The structure is chosen such that the 
gap changes from Ai to the lower value A 2 <C Ai 
when crossing junction 1, but it remains constant at 
A2 across junction 2. The gap has then also to change 
from A2 back to Ai in the loop part of the system, for 
example, via an extra large Josephson junction, which 



is irrelevant for the parity transitions in the small is- 
land. The interference part (8) is now generalized to 
Pi{2) = (1/4) [1 ~ cos< Pl{2) A m A 2 /E k E l + tk^/EkEi]. 
For parameters such that Ai + A 2 > SE > 2A 2 Cooper- 
pairs can be broken in a tunneling process across junction 
2. This leads to an extra term for the junction 2 tunnel- 
ing of the form 

-1 /"OO />oo 

f 2 = „ — 2 / / d&d&JMf*,6) 

JiT2<2 J — 00 J — 00 

x [1 - - fifa)]S(E k +E,- 6E), (10) 

where P2 = 1/2 — P 2 . Assuming that tunneling across 
junction 2 dominates we obtain 

5E 

r T «f 2 «- ^(l-ecos^), (11) 

where e equals -1 at the threshold SE = 2A and decreases 
weakly for larger SE. In this setup it is essential that i?T2 
is large enough such that the switching rate does not ex- 
ceed the qubit splitting Ej. For this to be true we need 
very asymmetric circuit Eji 3> i?j 2 . Under these condi- 
tions also 5E has no magnetic flux dependence, and the 
modulation in the transition rate originates completely 
from tunneling across junction 2. Numerical results for 
the rates are shown in Fig. 3b. 

Conclusion. We have re-analyzed in a new setup 
the phase dependence of quasiparticlc tunneling across a 
Josephson junction. A split Cooper-pair box embedded 
in a superconducting ring provides an optimal system 
to study the effect, as only quasiparticlc tunneling can 
create transitions between even and odd island charge 
parity, and the flux dependence of this rate occurs practi- 
cally due to the quasiparticlc-pair interference. The tran- 
sitions can be seen as jumps in the phase of a reflected 
weak microwave signal applied to the gate. The possible 
presence non-equilibrium quasiparticles makes the effect 
observable at very low temperatures, although (especially 
with differing gaps Ai 3> A 2 ) it is always possible to in- 
crease the thermal tunneling rate by heating the system. 
Transitions due to thermal and non-equilibrium quasi- 
particles result in a positive sign for e, whereas processes 
provided by breaking of a Cooper-pair in a negative one. 
A possible broadening of quasiparticle states would also 
lead to a negative sign for e. Therefore, the detection of 
the discussed phase dependence could also confirm non- 
equilibrium quasiparticles as the main source of quasi- 
particle poisoning in superconducting electronics. 
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